Impact of Subleading Corrections on Hadronic B Decays 
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We study the subleading corrections originating from the 3-parton [qqg) Fock states of final-state 
mesons in B decays. The corrections could give significant contributions to decays involving an 
LJ or T)^'^ in the final states. Our results indicate the similarity of loK and utt" rates, of order 
5 X 10"'', consistent with the recent measurements. We obtain a2{B — > J/ipK) ~ 0.27 + 0.05i, in 
good agreement with data. Without resorting to the unknown singlet annihilation effects, 3-parton 
Fock state contributions can enhance the branching ratios of Kr)' to the level above 50 x 10~®. 
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The rare B decays allow us to access the Kobayashi- 
Maskawa (KM) mixing angles and search for new physics. 
Much progress in the study of B decays 07 S 13 has been 
recently made in QCD-based approaches. In the pertur- 
bative QCD (pQCD) framework, the importance of the 
weak annihilation effects in _B ^ Kn decays was first 
emphasized by S], where the annihilation contributions 
are almost pure imaginary and therefore could lead to 
CP asymmetry predictions different from the QCD fac- 
torization (QCDF) results Nevertheless, the QCDF 
study showed that the annihilation effects may play only 
a minor role in the enhancement of tttt, t:K branching 
ratios (BRs) l^. A recent QCDF fit to Kn,!^ rates Q 
indicated that even if the annihilation contribution is ne- 
glected, one can still get quite good fitting results pro- 
vided that the strange quark mass is of order 80 MeV. 

The annihilation effects might be much more impor- 
tant for VP modes, where P and V denote pseudoscalar 
and vector mesons, respectively. It has been pointed out 
that in the absence of annihilation effects, the 4>K BRs 
are w 4 x 10~^ H which is too small compared to the 
data ~ 8 X 10~^ p, «^]- Recently Belle observed a large 
LuK- rate, (6.7li;f±0.6) x 10"^ and luK-/ujtt- ~ 1 [||. 
Sizable uK results are also reported in new BaBar mea- 
surements 6] with ujK~-° ~ ujTT^ ^ 5 X 10"®. It is hard 
to understand the large strength of ujK rates from the 

theoretical point of view. The ratio ujK /ojtt" reads 



UjK /uTT « \Vcb/Vub\^ ifK/Uf 



an 



agr^ + 2r2(a3 + 05 + ag/A) + fsfKba 



ai + ria2 



(1) 



where ri 



f^Ffyf^A^-, r2 = (F/^^/.)/(Af-/K), 
r^' = 2m'^ / [mb{ms + m„)] is the chirally enhanced fac- 
tor with nis^u being the current quark masses, and 63 = 
b^i^KjUj) is the annihilation contribution defined in 0. 
The uj'K~ rate depends weakly on the annihilation effects. 
Without annihilation, since 04 and a^r^ terms in the 

— . . . . — 

ivK amplitude have opposite signs, the ratio uiK /cjtt" 

should be very small. A possibility to explain the data is 
that the annihilation effects may give the dominant con- 
tribution to ujK modes as shown in the QCDF fit 8] for 
B PP together with some B — > VP modes. (However, 



including the contributions from annihilation effects, the 
pQCD results read Br(# ^ ujK°) < 2 x lO"® 0.) This 

result hints that, to account for the large ojK'^ rate, the 
annihilation contributions to the BRs of all B KV 
modes should be over 80%. If it will be true, it should be 
easy to observe, for instance, the following simple rela- 
tion p+X-'O : pO^-'O : t^i^-'O « 1 : [l/V^? : (l/^/2)^ 
the same as their annihilation ratios squared. Neverthe- 
less, if the global fit is extended to all measured B PV 
modes, a small Kuj rate ~ 2 x 10^® will be obtained ^3 
and a reliable best fit cannot be reached. The present 
QCD approach seems hard to offer a coherent picture in 
dealing with B VP modes. 

In this letter we take into account the subleading cor- 
rections arising from the 3-parton Fock states of final 
state mesons, as depicted in Fig.^ to QCDF decays am- 
plitudes. We find that it could give significant corrections 
to decays with w, or Ty''-* in the final states. A simple rule 
extended to B ^ PP, VP modes is obtained for the ef- 
fective coefficients af^ with the subleading corrections. 
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where i — 1, • • • , 5, and are the Wilson coefficients 
defined at the scale Hh = ~ 1.4 GeV with A^ 

the momentum of the emitted gluon as shown in Fig.^b) 
and 



/3 



V2 



Oi|S-)qqg = 0.12 (3) 



in the SU(3) limit. Here d = S7''(l - 75)m 1*7^(1 - 75)6, 
and oig is the averaged fraction of the tt~ momentum 
carried by the gluon. For the loK amplitudes, the term 
0-3 + 0-5, which is originally negligible, is replaced by 
03 + 05 + (c4 — C6)/3 and the latter gives the significantly 
constructive contribution to the rates. It can thus help 
understanding the reason for the similarity between Klu 
and ir^Lo. On the other hand, the subleading corrections 
can contribute significantly to the processes with 77^'^ in 
the final states, for which the term 03 — as always appears 
in the decay amplitudes and becomes 03 — + {a + cg) fs 
after taking into account the corrections. We also get 
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(a) 
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FIG. 1: The contributions of the qqg Fock states of the (a) ui 
and (b) 7r~ mesons to the B~ — > ljit~ amplitude. 



a^i^J/ipK) fa 0.27 + 0.05z which is weh consistent with 
the data. The result resolves the long-standing sign am- 
biguity of Re(a2). 

Let us study the subleading corrections originating 
from the 3-parton Fock states of final-state mesons. Tak- 
ing the cjTT" mode as an illustration, there are two dif- 
ferent types of diagrams shown in Fig. ^ In the follow- 
ing calculation, we adopt the conventions ^ da + 



0123 



-1, and use 



i5,TM^, Gap = (l/2)e„;3^,G'^^ e 
the Fock-Schwinger gauge to ensure the gauge-invariant 
nature of the results 



(4) 



For Fig. n^a) where the contributions come from the 
3-parton Fock states of the uj, because of the V — A 
structure of the weak interaction vertex, the relevant 3- 
parton light-cone distribution amplitudes (LCDA) up to 
the twist-4 level are given by 

{uiPuj, X)\u{0)j^gsGa/3ivx)u{0)\0) 



V2 
1 



Va e 



ip^XVCtg I ( UJ 



+ 7:r^P'^M''-' -Paa^/3)(*(a^) - \, (5) 



where Va = da a da „ da g (5(1 



ttg), with 



ag, a„, Ug being the fractions of the uj momentum carried 
by the u-quark, it-quark and gluon, respectively. Here $ 
and 5" are the twist-4 LCD As. Note that all the compo- 
nents of the coordinate x should be taken into account 
in the calculation before the coUinear approximation is 
applied. The exponential in Eq. Q before the coUinear 
approximation is actually e'^'^n'^'", where kg is the gluon' s 
momentum, and the resultant calculation can be eas- 
ily performed in the momentum space with substituting 
Xa —{i/v){d/dkg). The result of Fig. ^a) is found to 
be 



(WTT \Oi\B )Fig.l(a) = ft 



Due to G-parity, $ and ^ are antisymmetric in inter- 
changing and a„ for the w, so that Eq. © vanishes. 

In Fig. mb) , we consider the emitted gluon which be- 
comes a parton of the pion. We first take G^u{vx) ~ 
G^i/(0)e™'^s and then adopt the coUinear approxima- 
tion fcj = cfgPTT in the final stage of the calculation, where 
oig is the averaged fraction of the pion's momentum car- 
ried by the gluon. The calculation is straightforward and 
leads to 



t:-\Oi\B 



/Fig. 1(b) 
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Tr 



(vk^l -I- up^ 



2V2/. 



^pZ{n-\dri^gsGa^.b\B- 



(7) 



where the lj mesons's asymptotic leading-twist distri- 
bution amplitude 4>i^{u) = 6uu has been taken and 
u = I — u. We have two unknown parameters 
p"{TT^\d'-f^^'-f5gsGafib\B~) and ctg needed to be deter- 
mined. First, let us evaluate p^{Tr~\d'-f'^'y5gsGaiib\B~) . 
The matrix element can be calculated by considering the 
correlation function: 



Ilaip,p + q) = i 



d*xe 



tpx I 



-{q)\T{j^^{x)3Bmn 



1 



nib rrig ^ (P + (lY 
X {■K~\d-1^'-1zgsGapb\B- [q + p)) 



(8) 



where jzp = dgsGa^l^l^b, js = bi-f^u, the ellipses 
denote contributions from the higher resonance states, 
which can couple to the current js, and the transition 
matrix element can be parametrized as 



= p„/_(p2) + (p„ + 2q„)/+(p2). 



(9) 



In the deep Euclidean region of [p + q)^ , as depicted 
in Fig. 121 the correlation function can be perturbatively 
calculated in QCD and expressed in terms of 3-parton 
LCDAs of the pion. 



jjQCD 



du 



ml - {P + uqy JO 



dan 



X [-2{p ■ q)f3^cj)3^ + Umb{(t>\\ - 2(t>i:)], (10) 

where u = a^ + a^, and the 3-parton pion LCDAs are 
defined by [lllil 

{■K{q)\d{x)gsG^^ (ua::)(TQ;375'"(0) |0) 

= 'ifaAipiif^gi^a - q^gf,a) ~ qa{qt,gup - qugt^p)] 

Va(t)3^e"^''^°"'+''°'^\ (11) 
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FIG. 2: The diagrammatic illustration to the correlation func- 
tion, Eq. 0. 
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X / Va{(t>ii+(j)j^y 
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Here (f)^^^ is a twist-3 DA, and ^ii are all of twist-4 



FIG. 3: (a) Form factor f-{in^) plotted as a function of the 
Borel mass squared . (b) f-(rn^) — du dagAf_ with 
M^ = 10 GeV^ The volume in the plot is equal to f-{ml). 
Here u — Od + ag. 



f± prediction by adopting sq ~ 37 GeV^ and sa 
(9-20) GeV^. The result is depicted in Fig.OIa). The re- 
sulting value is f^{ml) = ±(0.057± 0.005) GeV^, where 
the uncertainty comes from the sum rule analysis. We 
then get 



1/ 

1 + ^1,02(7^9 - 3) 
+ cj2,o(2 - ^adOiu - Sag + Sa^) 



pZ{iT-\dn^9sG^,b\B-) 



-1.6 GeV*. (16) 



t^i,i(3adaQ - 2ag + Sa^ 



Mc^^) = 30(5"aj(l-ag)[- + 2e(l-2ag)] 
^\\{a^) = 



,1 



l20S^adauag[- + e{l — Sag)] . (13) 



Since the quark's momentum after emitting the gluon is 
roughly of order 7715/2 and the emitted gluon's momen- 
tum is Ax agp^ {ag will be discussed below), we set We therefore obtain 



Next, we determine the value of ctg. For illus- 
tration, we plot in Fig. |3tb) the amplitude Af_ of 
the f-i'rn'^) sum rule versus Ug and u{— a^, + ctg) 
by adopting M^=10 GeV^, where Af_ satisfies /_ = 

du Jp" dagAf^ , i.e. the volume in the plot is equal to 
/_(77i^). The resultant form factor is dominated by the 
region where u > 60%, ctg < 30%. The averaged frac- 
tion of the pion momentum carried by the gluon is then 
estimated to be cfg = (/^ du dagagAf_)/ ~ 0.23. 



/Fig.l(b) 



0.13 



the scale for the separation of the perturbative and non- 
perturbative parts at /i/i — ^ k^mBj^ ~ 1.4 GeV. The 
corresponding parameters at the scale [i^ read: /stt = 
0.0032 GeV^ wi,o = -2.63, = 9.62, wi.i = -1.05, 
52 = 0.19 GeV^', £ = 0.45 [13. To calculate /±, the 
contributions of higher resonances in Eq. © are approx- 
imated by 



ImngCD 

s - (p + q)- 



;ds, 



(14) 



where sq is the threshold of higher resonances. Equating 
Eqs. JSJ) and (I1U|) and making the Borel transformation: 
B[rn?g — {p + qY]^^ = exp(— 77i^/Af^), we obtain the 
light-cone sum rule 



du 



2m~gfB Jq J a 



dag e*- 



— '/>37r - /rr — (^11 - 24>±) 



(15) 



and/+(p2) ^ -/_(p2)^ where A = u-{ml-p'^)/{so-p'^). 
Using the above parameters for LCDAs, = (4.7 ± 
0.1) GeV, and fs = 180 MeV, we obtain the stable 



and /3 — 0.12 which gives the correction to 
as defined in Eq. jSJl. We list Oj without and 
with the subleading corrections in Table ^ where 
the approximation —^/2ip"{u!\u'-f^^gsGatib\B^)/AQ'^ ~ 
pZ{TT^\d'y^'-f5gsGaf,b\B-)/F^'' has been made. Note 
that oe and ag do not receive subleading corrections. 

In the following analysis, the LC sum rule form factors 
and nis ~ 80 MeV are used. We will instead use a smaller 
Aq'' = 0.28 which is preferred by the ujtt~ data. The ojtt~ 
mode is ideal for extracting Aq'' since its rate is insen- 
sitive to annihilation effects. We find that the spectator 
parameter Xh — ln-^(l + pne^'^") is consistent with 
zero in the analysis. tIic reason is that since the specta- 
tor interaction with a gluon exchange between the emit- 
ted meson and the recoiled pseudoscalar meson of twist-3 
LCDA $(j is end-point divergent in the coUinear expan- 
sion, the vertex of the gluon and spectator quark should 
be considered inside the pion wave function, i.e., for this 
situation the pion itself is at a 3-parton Fock state. An- 
nihilation effects have been emphasized in <j)K studies Q . 
We adopt the annihilation parameters pa — 0.9, 4>a — 
which give Br{B- 4>K-) ~ 8.5 x 10"^, consistent 
with the current data. Here the annihilation paramet_er 
of VP modes is defined as 



In 22^(1 + PA e"^-*) Q 
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TABLE I: Values for ai for charmless B decay processes without (first row) and with (second row) 3-parton Fock state 
contributions of final state mesons, where 03-10 are in units of fO~^ and the annihilation effects are not included. 



ai a2 0,3 a4 as ae a? as ag aio 

1.02 + 0.014* 0.10-0.08i 26 + 26i -328 - 91j 1.2 - 30i -487 - 72i 0.7 + 0.3i 4.5 + 0.6i -89 - O.li -5.9 + 7i 
0.974 + 0.0142 0.25 - 0.08i -55 + 26'i -291 - 91i 112 - 30i -487 - 72i -0.3 + 0.3i 4.5 + 0.6i -88 - O.li -18 + 7i 
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FIG. 4: (a) Dash, solid, and dot-dash for B 
ujn'^ ,u)K~ and uoK ; (b) solid, dots, dot-dash, and dash for 
K~ r)' r]' , K*~ r), and K'^rj. Brs are in units of 10^''. 
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FIG. 5: (a) Solid, dash, dot-dash, and dots for B 
K fP , K~ , K p~, and K~ (P\ (b) solid, dot-dash, dash, and 
dots for B K*~ -K^ ,K*'^ Ti^ , p~ r] and p'r]' . Brs are in units 
of 10"^ 



and its the imaginary part is neglected since the BRs 
are insensitive to it. In Fig. EJa) we plot the BRs of 
WTT" and ujK modes versus 7 (= argl^*^). The results for 
7 ~ (60—120)° are in good agreement with data. At 7 = 

90°, it gives wtt", wK", to be 5.5,4.5,4.3, respec- 
tively, in units of 10~^. Without the contributions from 
3-parton Fock states of mesons, (t)K~ ,ujtt~ ,ujK~ ,ujK 
wiU become 11,3.9,3.1,2.9 (in units of 10"^). The 3- 
parton Fock state effects give constructive contributions 
to WTT, LoK modes, but destructive one to the (j)K mode. 

The corrections from 3-parton Fock states of the kaon 
also give a definite answer to the longstanding problem 
for a2{J/ipK). In the earlier study, to account for the 
experimental value [02 1, the parameter pn has to be > 
1.5 01 . As emphasized in passing, without fine-tuning 
Ph, we calculate the amplitudes from the 3-parton Fock 
states of the kaon. With the same procedure as shown 
above, we obtain Z+C'tIj^^) ~ ±0.08, /a ~ 0.14 and thus 

a|L = af + ci{nh)f3 - 0.10 + 0.05i + ci(^,,)0.14 = 0.27 + 
0.05i, where a*^ is determined up to the twist-2 order 
and the SU(3) approximation for has been made. The 



result for 03 is well consistent with that extracted from 
data. This solves the long-standing sign ambiguity of 
o.2{J /ipK) which turns out to be positive for its real part. 
Note that if Rc(a2) were negative, /s would have to be 
~ —0.3 which in turns would lead to (t)K ~ 20 x 10~^ 
and w7r-,wK~ 1 x lO^^! 

The subleading corrections could give significant con- 
tributions to the decays with r/'^ in the final states be- 
cause these decay amplitudes always contain the singlet 
factor 03 — 05. We plot the BRs of Krj',K*r} modes 
versus 7 in Fig. Efb), where « 0, the annihila- 

tion parameter for PP modes, has been used as it could 
give good fit results for K-k^tttt rates 0- We do not 
consider the singlet annihilation correction |l5j | in K-q' 
modes because it is still hard to determine at present. 
With (without) the subleading corrections, we see that 
K-i > i?V « 55(35), and K*-r] > K*"r] « 24(20), 
in units of 10~^. The corrections give 70% and 25% en- 
hancements to Krj' and K*ri rates, respectively. Note 
that with (without) the corrections, K*ri' < 1 x 10^^ 
(> 4 X 10"^) and Kr]<lx IQ-^ (> 1 x 10"^). In pQCD 

calculation 16j, it seems that K ri'{— 41 x 10~^) was un- 
derestimated while K'^ri{= 7 x 10^^) overestimated com- 
pared to the data0,l3. Within the QCDF framework, 
by only considering the two-parton LCDAs of mesons, 
Beneke and Neubert obtained K*rj ~ 13 x 10^^, just 
half of the experimental value, but with a huge error. 

For further comparison with other calculations, we plot 
Kp, K*!^, p~ri'''^ modes in Fig. |S1 The subleading con- 
tributions to these BRs are < 15%. At 7 = 90°, we 
ha.yeK°p°,K-p+,Tfp~-,K-p° = 6, 7, 7, 2 (xlO^^), and 
K*~TT+,K*°TT°,p-r],p-r]' = 9,2,6,4 (xlO^^). In 0, 
to fit (pK, LoK rates, the annihilation effects dominate 
the decay amplitudes and the form factors , A^'^ 
are rather small, such that it will lead to p^K^-^ : 
p^K'-^ : LuK-'° « 1 : (l/^/2)2 : (1/72)^, while the 

pQCD results for rates are k'^ p'^, K^p^, k'^ p~ , K~ p^ — 
2.5,5.4,3.0,2.2(xl0-6) 0. 

In conclusion, we have calculated the contributions 
arising from 3-parton Fock states of mesons in B de- 
cays. We find that the contributions could give signif- 
icant corrections to decays with an lo or Ty*-'-* in the fi- 
nal states. Our main results for 7 « (60 — 110)° are 

ijjTT^ ^ljjK~ ^LjjK'^ ~ 6.0, (6 5), 5.1, respectively, in units 
of 10^^, while the previous QCDF global fit to VP modes 
and the pQCD results gave smaller loK BRs of order 

< 3 X 10-6 H US- We predict that i?°pO - ujK and 
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K /K~p^ 3. Including the corrections we obtain 
0'2{J/tpK) « 0.27 + Q.Obi which is well consistent with 
the data. The sign of Re(a2) turns out to be positive. 
Without resorting to the unknown singlet annihilation 
effects, 3-parton Fock state contributions can enhance 



Kri' to the level above 5 x 10 ^. 
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